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We develop a renormalized perturbation theory for the dynamics of interacting Brownian particles, 
which preserves the fluctuation-dissipation relation order by order. We then show that the resulting 
one- loop theory gives a closed equation for the density correlation function, which is identical with 
that in the standard mode coupling theory. 
PACS 64.70.Pf 

As a first-principle approach, the mode coupling theory (MCT) [2, HOI has not only enjoyed considerable suc- 
cess in explaining the slowing down of the weakly supercooled liquids, but also had enormous impact on the area 
by stimulating further experiments, simulations, and other theoretical developments. However, the foundation 
of the theory needs much to be desired since the theory is fraught with uncontrolled approximations. It is helpful 
if one can develop systematic field-theoretical treatment for model system containing smallness parameter. The 
earlier field-theoretic formulations [J, [5[ of MCT are found to be incompatible with the fluctuation-dissipation 
relation (FDR) It was thus an urgent task to develop a consistent field-theoretic formulation capable 

of describing the dynamics of glass-forming liquids, for which a systematic perturbation expansion preserving 
the FDR is possible. Quite recently, Andreanov, Biroli, and Lefevre (ABL)Q provided a remarkable insight 
into this problem by focusing on the symmetry properties of the action integral under time-reversal (TR). In 
particular, ABL has identified the transformation of fields under TR, which leaves the action invariant. The 
FDR naturally follows from the TR transformation combined with causality. Moreover, the nonlinear nature 
of the TR transformation is shown to be the underlying reason why the renormalized perturbation theory (i.e., 
the loop-expansion) does not preserve the FDR order by order. By introducing a new set of auxiliary fields 
to linearize the time-reversal transformation, ABL have attempted to develop a FDR-preserving field theory. 
Although ABL's work is a remarkable step forward, the one-loop results of ABL's theory are found to have some 
pathological features. The equation for the nonergodicity parameter gives nontrivial results even for noninter- 
acting Brownian systems, which is suspicious and should be examined carefully. Furthermore, their vertex is 
ill-behaved, leading to the divergence of the memory integral. We tend to believe that these ill-behaved results 
are intimately connected to their linearization scheme. 

In this Communication, by proposing a simpler but crucial linearization scheme, we show that the one-loop 
result in the FDR-preserving renormalized perturbation theory yields a closed dynamic equation for the density 
correlation function, and show that this closed equation turns out to be the same as that in the standard MCT. 
We thus have established the precise relationship of the FDR-preserving field theory with the standard MCT. 

We start with the following Langevin equation for the density field p(r, t) of interacting Brownian particles 

Mr, *) = V • (p(r, *)V^Jy) + V(r, t) (d) 

where the Gaussian thermal noise rj(r,t) has zero mean and variance of the form 

< r)(r, t)r](r', t') >= 2TV • V (p(r, t)S{r - r')S(t - t')) (.2) 

Note that the noise correlation depends on the density variable, i.e., the noise is multiplicative. In (Q]), F[p] is 
the free energy density functional which takes the following form: 



F[p]=T J drp(r)(ln^-l)+i J dr J dr' Sp{r) U(r - r') 6p{r') 



(.3) 



where 5p(r, t) = p(r, t) — po is the density fluctuation around the equilibrium density po. In (Q||) the first term is 
the ideal gas part of the free energy, F i( i [p] , and the second term the interaction part of the free energy, F int [p] . 
Using Ito calculus, Dean Q has derived the above nonlinear Langevin equation for the density field of system of 
interacting Brownian particles with pair potential U(r). Earlier, Kawasaki [9( has obtained the same Langevin 
equation with U(r) replaced by — Tc(r), c(r) being the direct correlation function, by adiabatically eliminating 
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the momentum field in the fluctuating nonlinear hydrodynamic equations [J] of the glass-forming liquids. For 
this case, the free energy density functional (Q|]) takes the Ramakrishnan-Yussouff form. 

We consider the corresponding action integral S[p, p] which governs the stochastic dynamics of the coarse 
grained density variable 

S[p,p}= J drJdt{p[-d t p + V-[pV S -^)]+Tp(Vp) 2 } (.4) 

where the field p is pure imaginary, and the last term comes from the multiplicative thermal noise. (A similar 
action was given in [5| but with real p.) It is a crucial observation of ABL to recognize that the above action is 
invariant under the TR field transformation: 

P(r,-t) = p(r,t) 

p(r,-t) = -p(r,t) + ^ I ^ J) (.5) 

The FDR follows from the above TR transformation. The response function R(r,t;r't') is defined as a link 
between induced density change A < p(r,t) > and an external infinitesimal field h e (r',t'): 

A < p(r, t) >= J dr' J dt' R(r, t; r't') h e (r', t') (.6) 

where <■••>= J Vp J Vp (• • ■) exp(<S[p, p]). To find the induced change of density, one should add the contri- 
bution of the external field to the free energy F, AF = — f dr f dtSp(r, t)h e (r, t). It is straightforward to show, 
by obtaining the induced density change, that the response function R(r, t; r't') is given by 

R(r, t; r't') = -(p(r, t)V ■ (p(r', t')V'p(r', t')) ) (.7) 

Note that the response function is not given by the conventional response function — < p(r, i))V' 2 /5(r', t') > 
which is the true response function for the Langevin equation with additive noise. Instead, the response function 
is given by (JT7J) due to the multiplicative nature of the Langevin equation noise for the density fluctuation [|| 0] • 

Using the identity (p(r,t)SS/Sp(r ! ,t') \ = and the TR transformation |3]), one obtains the FDR 



- ^d t G pp (r -r',t- t') = -R(r - r' , t' - t) + R(r -r',t- t') (.8) 

where G pp (r -r',t- t') = (Sp(r, t)Sp(r', t')). Since R(r - r', t' — t) = for t > t' due to causality, @ gives the 
standard form of the FDR 

R(r -r',t- t') = -B(t - t')±d t G pp (r -r',t- t') (.9) 

where Q(t) is the Heaviside step function. 

Using the form of the free energy given in ([75]) . one can explicitly write 5F/Sp as 



T 5p{r, t) po po 

1 SF int [p] 
TS P (r,t) T 



l - J dr'U{r-r')5p(r',t) (.10) 



where f[Sp(r,t)] is the contribution of the non-Gaussian part of Fid[p], and is given by an infinite series 
f[Sp(r, t)] = — J2^=2 n(~ ^/°( r ' / Po) ■ Using (j.lOp one can explicitly write the TR transformation (JT5J) as 

P(r,-t) = p(r,t) 

p(r,-t) = -p(r,t) + f[5p(r,t)]+K*Sp(r,t) (.11) 

where K* is convolution with the kernel K(r) = (6(r)/po + U(r)/T). Note that (j.lip is nonlinear due to the 
non-Gaussian contribution f[5p(r,t)] to the ideal-gas part of the free energy [l(|- ABL has shown that this 
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nonlinear nature of the transformation is the origin of the incompatibility of the renormalized perturbation 
theory with FDR. Naturally, this incompatibility would be resolved if the transformation (jTTTJ) is made linear 
by ignoring f[Sp(r,t)]. However, as shown below, in this Gaussianizcd case, the nonlinear term generated by 
dynamics (the second term in (|.13[) below) would give rise to the spurious contribution to the one- loop result, 
incorrectly yielding a nontrivial result [f| @, HH even in the noninteracting system. 

As the most natural way to make the transformation (j.lip linear, we introduce a new field 9(r,t) so as to 
satisfy the nonlinear constraint 

00 i x 

9(r,t)=f[Sp(r,t)] = -J2-(-^) (-12) 



Note that our approach here differs from that of ABL in that whereas ABL defines the new variable as the 
functional derivative of the full free energy with respect to density: Q abl (r, t) = SF/Sp(r,t), we limit the new 
variable 9(r, t) only to the nonlinear part of it. Using the constraint (|.12D . we obtain the ideal-gas contribution 
to the body force as 



V • (pV— —) = TV 2 p + — V • (5pVp) + p TV 2 9 + TV • (5pV9) (.13) 

V OP ) Pn 



where the first (the last) two terms are the contributions from the Gaussian (non-Gaussian) parts of the ideal- 
gas free energy. Since due to cancellation of the two nonlinear effects, the entire ideal-gas contribution to the 

dynamics is pure diffusion, i.e., V • (^p\7 S Fid / S pj = TV 2 p, the sum of the last three terms in (|.13|) should vanish 

if the constraint (|.12[) is taken into account: 

—v • (SpS7p) + p Tv 2 e + tv ■ (Spve) = (.14) 

Po 

As shown below, this nonperturbative cancellation plays a crucial role in obtaining a closed equation for the 
density correlation function alone. 

Incorporating the new variable 9(r, t) and its conjugate 6(r, i), the action can now be explicitly rewritten 

as 

sty] = s a [i/j] + s ng [i>] 

S g [i>] = J dv J dt {p [ - dtp + TV 2 p + p TV 2 9 + p V 2 U * Sp\ +Tp (Vp) 2 + <96»} 
S ng [ip] ee J dr J dt{p V • (dpVU * <Sp) + — V • {5pVp) + TV • (6pV6) 

+ T6p(\7p) 2 -9f(6 P )} (.15) 

where tp denotes the entire set of the fields collectively, and the full action S[i[>] is separated into its Gaussian 
part S g [tjj] and non-Gaussian part S ng [i/j] . Now the actions S g [ip] and S ng [tp] are separately invariant under 
the following linear TR transformation 





-t) = 


P(r,t) 


/5(r, 


-t) = 


-p(T,t) + 9(r,t) + K*Sp(T,t) 


9(r, 


-t) = 


9(r,t) 


6(r, 


-t) - 


9(r,t)-d t p(r,t) 



(.16) 

It is easy to show that the modulus of the associated transformation matrix is unity. Though the three underlined 
terms in (|.15p vanish when summed together, their presence is crucial for the actions S g [tp] and S ng [tp] to 
be separately time-reversal invariant. Note also that the separate invariance of the actions under the linear 
transformation (|.16p is not tied to the form of the constraint (|.12[) . This separate invariance of <S g [^>] and S ng [ip] 
enables us to construct the FDR-preserving renormalized perturbation theory from these actions. 
With the new action (|.15[) , it is easy to show that the response function is given by 

R(T,t;r',t') = ^5p(r,t)§(r',t')) (.17) 
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One can then obtain the FDR fcEfy by taking correlation of the last member of (|.16l) with Sp(r,t)/T. 

Incorporating the above new set of variables into the theory, we are now ready to develop a renormalized 
perturbation theory which preserves the FDR order by order. We first consider the noninteracting case (U = 0) 
and show that the noninteracting part of the action, Sid[tp] = S[ip; U = 0], yields the dynamic behavior consistent 
for the noninteracting system. We begin with the following identities 

where 1 = (r, t) and 2 = (0, 0). The first identity can be written explicitly as 
= (^( 2 )^#) = ( - | + TV 2 )G PP {1 2) - 21> V 2 ( / 5(l),5p(2)) - 2T(^(2)V • (^(l)Vp(l))) (.19) 



where we used the fact that the sum of the three underlined terms in (|.15j) vanishes. Similarly, using the second 
identity in (|.18[) . we obtain 

= (*P(2)^|#) = PoTV 2 {p{l)5 P {2)) + (9(l)Sp(2)) + T(Sp(2)V ■ (^(1)V/5(1))) (.20) 

In (p20f , (p(l)Sp(2)) = aa(6(l)Sp(2)) = for t > due to the causality, and hence we obtain 

(Sp(2)V- (fy(l)Vp(l))) = for <>0 (.21) 
Using (l.2ip and causality, we obtain from (|.19[) 

^G pp (r,i) =TV 2 G w (r,t) for i > (.22) 

This result is expected for the non-interacting system. 

There are five nonlinear terms in the full action (|.15p where the two underlined nonlinear terms (T/po)pV • 

(SpVp^ and TpV • (jpVtf) are shown to cancel the linear diffusion term ppoTV 2 9. Then the remaining 
nonlinearities are pV • \ 5pVU * 5p^j, T5p(V pf , and 0/[<5p]. The first two come from the particle interaction 

and the multiplicative thermal noise, respectively. The last cubic nonlinear term ^8(6 p/ po) 2 , the only one 
contributing in the one loop order, comes from the non-Gaussian part of the ideal-gas free energy. In order to 
analyze the effect of these three nonlinear terms, one should examine the structures of the relevant self-energies. 
Let us write 

SM = ^ T (l)-G^(12)-m 

SnM = ^(123)^(l)VX2)^(3) + 0(l)f;i(-^) ri (.23) 
3! po J 

where and i/j t (1) are respectively column and row vectors of the four fields p, p, 9, and 9, and the term 
with n = 2 in the summation is incorporated into V^123). The unperturbed inverse matrix propagator G ~ 1 (12), 
can be read off from the action (|.15[) as 



Go 1 (12) = 



( Dx5{12) + p V 2 [/(12) 0\ 

D x 5{12) + p V 2 U(12) -2T Po V\8{12) 7> V 2 <5(12) 

7> V 2 <5(12) J(12) 

V (5(12) / 



(.24) 



where D x = {-d/dt x + TV 2 ) and D 1 = {d/dt l + TV 2 ). V(123) in fact is a collection of vertices V aiaaaa (123) 
with the a's standing for four field types, 1, 2, 3 ■ • • for space-time coordinates, and "0(1) etc stand for ip ai (1) etc. 
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Repeated thin (thick) numbers indicate space-time integrations and summations over the field types (space-time 
integrations only). We give expressions for the vertices V r ctlQ2Q3 (123): 



V$L(123) = -— Vii(12) • V 3 <5(23) + Vi<J(13) • V 2 5(23) 



T 



Po 

V%£(123) = ~[Vi«(12). V 3 ^(23)+Vi5(13)- V 2 #(23)]<5(i 2 -t 3 ), 
U ppe (123) = -TVi • V 3 *(123), V P pp{123)=TV 2 S(12)-V 3 6(13), 
V § PP ( 123 ) = ^^(12)5(23) 



.25) 



where we have separated the vertex Vpp P (123) into the ideal-gas contribution V~p„(123) and the interaction 
contribution VI™* (123). In (|.24j) and (|.25[) the space-time delta function is defined as 5(12) = S(ri-r 2 )5(ti~t 2 ). 

The dynamic equations for the correlation and response functions are formally given by the matrix Schwinger- 
Dyson (SD) equation 

G : (13) • G(32) - £(13) • G(32) = 5(12) (.26) 

In order to obtain the equation of motion for the density correlation function, we compute pp-element of the 
SD equation by taking 1 = (r, t), 2 = (0,0), and 3 = (r s ,s). It is now convenient to introduce the Fourier 
transform in space 

£(r,i) = / S(k,t)e 4k - r 



(.27) 

etc. where J k = J (ik/(27r) 3 . Using (|.24| and causality of the response functions, we obtain for t > 

[Gq 1 ■ G] pp (k,t) = -(d t + p Tk 2 K(k^Gp P (k,t) - Po Tk 2 Gg p (k,t) (.28) 
Likewise using the causalities of the self-energy functions and the response functions, we obtain for t > 

[E-GL (k,t) = 



ds 



ds 



Sp P (k, t - s)G pp (k, s) + Epe(k, t - s)Ge P {k, s) 



Spp(k, t - s)G pp (k, s) + S p g(k, t - s)G §p (k, s 



.29) 



where the upper limits of the time integration are due to the causality of the self-energy functions in the first 
two integrations, and to the causality of the response functions in the last two integrations. It can be shown 
that (1.291) can be simplified using the TR properties of the self-energies and the response functions as 

ft 



[E.GL(k,t)= - 



S^(k, t - s)d s Gp P (k, s ) + if (k)£ pp (k, t - s)G pp (k, s) 



/ ds S pp (k,t - s)G ep (k, s) 
Jo 



.30) 



Using (|.28p and (|.30|1 , the pp-element of the SD equation for t > can now be explicitly written as 



d t G pp (k,t) 



p a Tk 2 K(k)G pp (k,t) - p a Tk 2 Ge p (k,t) 



ds 



E p§ (k, t - s)d s G pp (k, s) + K(k)E P p(k, t - s)G P p(k, s) 

/ ds Spp(k, t - s)G 9p (k, s) 
Jo 



.31) 



One can likewise obtain ^p-element of the SD equation as 



G flp (k,t) = £^(k,0)G w (k,t)- f dsX §§ (k,t- s )d s G P p(k, 

Jo 



s) 



ds 



K(k)E § Jk, t - s)Gpp(k, a) + E fo (k, * - s)Gg p (k, a) 



.32) 



We now take the two crucial steps which can give a closed equation for G pp (v,t) in one-loop order: 



G 



The equations (| .31 [) and (|.32[) imply that the integral involving G$ p (r,t) in (| . 3 1 1) are of higher order and 
hence can be discarded in the one- loop order (note that [Go}e p — from (l.24[) ~). 

As we alluded earlier, using the cancellation of the three underlined terms in (|.15p . one can eliminate 
— poTk 2 Ge P (k, t) in ().31|) . together with the two vertices Vi p and V PP e appearing in (|.30|) . It should 
be emphasized that the cancellation is a nonperturbative effect which is required to hold only when the 
constraint (|.12p is employed. Therefore, we discard — p^Tk 2 Gg p (k, t) and retain only those terms arising 

7i%, V ppp , and V §pp 



from the remaining three vertices Vg"*, V ppp , and Vx in evaluating the one-loop self-energies. 



Taking these steps we obtain the equation of motion for the density correlation function up to the one-loop 
order as 



8 t G pp (k,t) = - Po Tk z K(k)G pp {k,t) 

+ [ ds \E &§ (k,t-s)d s G pp (k,s)+K(k)i: ?p (k,t-s)G pp (k,s)\ (.33) 
Jo 1 J 

In the absence of particle interaction (U = 0), the equation (|.33p reduces to the diffusion equation (|.22p since, 
as shown below, the self-energies in (|.33[) involves the vertex V~" p which contains the interaction potential U. 

Now we compute the one-loop self-energies E p ^(k, t) and E pp (k, t). With the three surviving vertices Vi pp , 
V ppp , and Vg , we find that there is the only one nonvanishing diagram for £.g(12), which is given by 

£^(12) = Vif p (134)V 8pp (256)G pp (35)G pp (46) (.34) 
Likewise, there are three nonvanishing contributions to the self-energy £ pp (12): 

E pp (12) = EW(12) + sg(12) 

Eg (12) = V^p(134)V^p(256)Gpp(35)Gpp(46) 

Eg (12) = ^; p 4 (134)V/ ppp (562)G pp (35)G pp (46) (.35) 

Eg in (j.35[) comes solely from the interaction contribution to the body force, and Eg from the cross- 
contribution of the interaction and multiplicative thermal noise. By multiplying the second equation in the 
transformation (|.16p by Sp(0, 0), taking average and using causality, one obtains 

G pp (r, t) = Q(t)(K * G pp (r, t) + G p9 (r, t)) (.36) 

When (|.36p is substituted into (|.35p . the correlation function G p g(r, t) will make null contribution in the one- loop 
order. Using this fact one can rewrite (|.35|) as 

Eg (12) = ^(134)y-;(256)G pp (35)G pp (46) 

Eg (12) = ^;*(134)^(562)G w (35)ir*G w (46)e(t 4 -t 6 ) (-37) 

Therefore the equations (| ,33[) - (| ,34[) and (|.37[) consist of a closed equation for the density correlation function 
G pp (r,t) alone. 

The Fourier components of the self-energies arc now computed as 



E (k,t) = --L f l/(k,q)G pp (q,i)G pp (k-q,t) 
Eg(k,t) = / F 2 (k,q)G pp (q,t)G pp (k-q,t) 

Jq 

Eg(k,t) = -l^y 2 (k,q)G pp (q,t)G pp (k- q ,t) 



T 
2p~ 



-k 2 / T/(k, q )G pp (q,0G pp (k-q,t) (.38) 

'0 Jo 
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where V(k, q) = [(k • q)t/(q) + k • (k — q)U(k — q)] . Arranging the sum of the self-energies Yr?~ with i = 1, 2, 
one can rewrite Epp(k, i) as 



E w (k,t) = k 2 E P p(k,i) 



E P p(k,i) 



V 2 (k,q) 



T 
Po 



-V(k,q) 



G P p(q, t)G pp (k - q, i) 



(.39) 



with V(k,q) = V(k,q)/|k|. 

In order to appreciate what the closed equation (|.33p implies, let us for the moment ignore the contribution 
of the self-energy E^(k, i) in (|.28p . Then (33) would become 



d t G P p(k, t) = - Po Tk 2 K(k)G pp (k, t) + k 2 K{k) / ds S P p(k, t - s)G pp {k, s) 



(.40) 



One can see from (j.40[) that the effect of the nonlinear contribution Epp(k, t) is to renormalize the 'bare' life 
time T (k) = (p Tk 2 K(k)) 1 || as follows. Defining the Laplace transform Gp p (k, z) = J °° dt e~ zt G pp (k, t), 
etc., we obtain the Laplace transform of (|.40p as 

(.41) 



G^ p (k, z) = G pp (k,0)- z + T^(k,z) 



where wc defined the rcnormalizcd life time as 

1 ,-. 



)(k)(l 



th(M) = r (k)( 1 - ^^p(k,z) ) = r (k)(l 



the last equality holding in the one- loop theory. Substituting (|.42[) into (|.4ip , one obtains 



G^ p (k, z)=G P p (k,0) 
The corresponding equation in time domain is given by 



z + 



^^(k) 



l + ^ p (k,z)/ Po T 



d t G pp (k, t) = -p Tk 2 K(k)G P p{k, t)+ J^ds E P p(k, t - s) ( - a s G P p(k, s)/p T 



(.42) 



(.43) 



(.44) 



The equation (|.44p is in fact obtained by replacing G pp (k, s) in the convolution integral in (| .40[) by 
( — <9 s G pp (k, s)/po7 1 k 2 i^(k)), which is valid up to the one-loop order. 

It is clear that the so far ignored self-energy T,~g(k, t) makes a new contribution to the renormalization of the 
bare life time ro(k), which comes from the consistency requirement of the perturbation theory with the FDR. 
Now restoring the contribution of Y,~q (k, t) and adding it to (1.44|) , we obtain 

d t G P p(k,t) = -p Q Tk 2 K(k)G P p(k,t)- f ds-Z MC (k,t-s)d s G pp (k,s), 

Jo 



->MC 



(M) = (-s^(k,t) + — v,.. Ik ./i 



l 

PoT' 



2 Po T 



V 2 (k,q)Gpp(q,t)G pp (k~q,t) 



.45) 



This equation reduces to the closed dynamic equation for the density correlation function in the standard MCT 
if U(k) is replaced by — Tc(k). The self-energy E-g(k, t) and the part of the self-energy EL ; (k,t) which is 
multiplied by poT cancel against each other to yield the last line in (|.45p . 

The equation for the non-ergodicity parameter (NEP) f(k) defined as /(k) = G pp (k,t — > oo)/G pp (k, 0), can 
be obtained from (1. 451) as 



-<MC 



(k,oo) 



1 



/(k) 

1 - /(k) p Tk 2 A'(k) 2 Po T 2 k 2 



V 2 (k, q)5(k)5(q)5(k - q)/(q)/(k - q) 



.46) 
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In obtaining (|.46p . we used the inverse relationship between K(k) and the static structure factor 5(k): poK(k) = 
(l + p U(k)/T) = (l - Poc(k)) = S ,_1 (k). The NEP equation (T46| is similar to the corresponding one in ABL, 
except that there is no memory kernel involved in ABL. It naturally occurs from (|.46[) that /(k) = is the only 
solution for the noninteracting case (U = 0). But this feature is absent in ABL's theory. 

In summary, in order to obtain the FDR-preserving perturbation theory, we here proposed a simpler but 
crucial linearization scheme in which the new variable is employed to take care of only the nonlinear part of the 
TR transformation. We then recognized that there is a characteristic nonperturbative cancellation effect in the 
theory. This feature enables us to obtain in the one-loop theory a closed dynamic equation for the two-point 
density correlation function. This equation is shown to be the same as that of the standard MCT. 

Having established the relation of our field-theoretical treatment with the standard MCT, we are now at the 
starting point to venture into ambitious tasks such as higher-order loop calculations and multibody correlations. 
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